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Abstract: We consider a test for the hypothesis that the within-treatment 
variance component in a one-way random effects model is null. This test is 
based on a decomposition of a [/-statistic. Its asymptotic null distribution 
is derived under the mild regularity condition that the second moment of the 
random effects and the fourth moment of the within-treatment errors are finite. 
Under the additional assumption that the fourth moment of the random effect 
is finite, we also derive the distribution of the proposed [/-test statistic under a 
sequence of local alternative hypotheses. We report the results of a simulation 
study conducted to compare the performance of the [/-test with that of the 
usual F-test. The main conclusions of the simulation study are that (i) under 
normality or under moderate degrees of imbalance in the design, the _F-test 
behaves well when compared to the f/-test, and (ii) when the distribution of 
the random effects and within-treatment errors are nonnormal, the [/-test is 
preferable even when the number of treatments is small. 

1. Introduction 

Consider the one-way random effects model 

(f .1) Y,.i = ^ + 6, + ey , i = f , . . . fc, j f , . . . , (> 2) 

with bi and e^j, i = l,...,k, j = l,...,ni, denoting independent random vari- 
ables with null means and variances and a^, respectively. Thus, /i is the mean 
response, bi represents the random effect of the i-th treatment, represents a 
random measurement error associated with the j-th observation obtained under 
the i-th treatment, and and are between- and within-treatment variance 
components, respectively. 

In general, data analysis based on such a model focuses on the estimation of 
and on testing the one-sided hypothesis (of no treatment effects) 

(1.2) Ho : cTfe = vs Hi : > 0. 

Inference about variance components in random effects models, and more generally 
in linear mixed models, has a long history in the statistical literature. In this con- 
text, Searle et al. [21] provides an excellent overview of estimation and prediction 

* Supported by Conselho Nacional de Desenvolvimento Cientifico c Tecnologico (CNPq) and 
Fundagao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP), Brazil. 

^ Dcpartamcnto de Estatistica e Matematica Aplicada, Universidade Federal do Ceara, Brazil, 
e-mail: juvencioSuf c.br 

^ Departamento de Estatistica, Universidade de Sao Paulo, Brazil, e-mail: jmsingeraime.usp.br 

■^Department of Econometrics and Business Statistics, Monash University, P.O. Box 197, 
Caulfield East, Australia 3145, e-mail: Mervyn.Silvapulle@buseco.monash.edu.au 

AMS 2000 subject classifications: Primary 62F03; secondary 62F05. 

Keywords and phrases: martingales, one-sided hypotheses, one-way random effects model, re- 
peated measures, [/-statistics, variance components. 



197 



198 



J. S, Nobre, J. M. Singer and M. J. Silvapulle 



while Khuri et al. [l.!] and Demidenko [6] present an extensive review of hypothesis 
testing; Silvapulle and Sen [M>] consider a comprehensive treatment of testing under 
inequality constraints. 

Under the additional assumption that bi and Cij follow normal distributions, the 
usual F-statistic for testing (1.2) is 

(1.3) SQ(6)/(fc-l) 



are, 



SQ{e)/{n-k)' 

where SQ(6) = E ti n^{Y,. - Y.^ and SQ(e) = E ti - E ti n,F^ 

respectively, the between- and within-trcatment sums of squares. The i^-statistic 
(1.3) follows a central F distribution with fc — 1 and n — k degrees of freedom when 
Tio : cr^ = is true. In the balanced case (i.e., when ni = ■ ■ ■ = Uk) the test is 
uniformly most powerful invariant (UMPI) . This optimality property does not hold 
in the unbalanced case. Details may be found in Khuri et al. [13], for example. 

We propose an alternative test based on the decomposition of [/-statistics con- 
sidered by Pinheiro et al. [19] in a nonparametric setup. Although it is not an exact 
test, it has good properties for moderate sample sizes and does not require the 
normality assumption. It is also valid for hcteroskcdastic random effects. 

The class of U -statistics has its genesis in the papers of Halmos [1 1 ] and Hoeffding 
[1 2] and is well known for its simple structure and for the weak assumptions required 
for its use in statistical inference. It also provides a unified paradigm in the field of 
nonparametric Statistics and has been used in many applications, as illustrated in 
Sen [23, 24, 25, 26, 27, 28, 29, 30, 31] and Sen and Ghosh [32], among others. The 
related theory is available in many sources, among which we mention Scrfling [M], 
Sen and Ghosh [32], Lee [14] or Sen and Singer [33]. 

The test is derived under the assumption that E[e^j] < oo and thus, accommo- 
dates a large class of distributions (not necessarily continuous) as the sources of 
variation in model (1.1). 

In Section 2, using a martingale property, we describe the decomposition of the 
[/-statistic that provides the key to the proposed test. In Section 3, we present 
simulation studies designed to evaluate the efficiency of the proposed test in differ- 
ent situations and to compare its performance with that of the exact i^-test. We 
conclude with a brief discussion in Section 4. 



2. Testing for variance components 

Consider the function g{x,y) = {x — 2/)^/2 and note that, under model (1.1), 

E[5(y,„y,,0] = E[(e,, - e,,0']/2 = crl 

Therefore, an unbiased estimator of (T^, based only on the Ui observations obtained 
under the i-th treatment is given by the following [/-statistic 

71 i 

(2.1) = (n, -1)-1E(^^.7- l=l,...,k. 

i=i 

Since E[{b^ - b^,){e,j - e,,f)] = 0, it follows that ¥.[g{Y,j ,Y,, j,)] = {2al + 2al]/2 = 
cr^ +(Tg . Therefore, an unbiased estimator of crj^ + <Tg , based only on the observations 
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obtained under treatments i and i' , is given by the following generalized [/-statistic 
of order (1,1) 

(Y- ■ - Y, ■,'|2 

(2.2) Uu' = [n,n,')-^Y.ll ' 2 ' 1 A:. 

i=i j'=i 

Letting n ~ z^^^i the lexicographically ordered observations, 

^11 , ■ • ■ J Ylnx , ^21 , ■ • ■ , ^2n2 : ■ • ■ i ^fcl j ■ • ■ i ^fcrifc ; 

may be re-expressed as 

Yl, . . . , Yn-^ , • • • , Yni+n2 i ■ ■ ■ i -^iiH hrifc-i-l-li ■ ■ ■ ; ^ni 

where the first ni observations relate to treatment 1, the next n2, to treatment 2 
and so on. The uniformly minimum variance unbiased estimator (UMVUE) of the 
variance of the observations is given by the [/-statistic 



2' 

1 < r < s < n 



which can be rewritten as 



1=1 ^ ' l<i<i'<k 



i=l ^ ' l<i<i'<k ^ ' 

highlighting its nature of a linear combination of generalized [/-statistics. The first 
and second terms in (2.3) correspond, respectively, to the within and between- 
treatments components. Since 

niijii - \)l{n(n - 1)} = {rii/ji} - {ni(n - ni)\l{n(n - 1)}, 

and 'Ylii'^i "-i' ~ n- ~ "n-i: we can rewrite the first term in (2.3) as 

E'^ ni{ni - 1) _ jT ~ jj 

n{n-l) ' ~ ^ n'^ ^ nin^X)^' 

i=l ^ ' i=l i=l ^ ' 

= t-u^- E -i^.u.- j: ^^[/., 

^ — ' n ^ — ' n ri — 1 ^ — ' nm— 1 

i=l l<i<i'<k ^ ' l<i<i'<k ^ ' 

leading to the decomposition 

(2.4) Ul = VK„ + B„ 

where 

h 

Wr, = y2^U, and B„= V {2[/,,, -[/,-[/,} . 

^ — ' n ^ — ' nln — 1 

4=1 l<'i<i'<fc ^ ' 
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This includes the classical ANOVA decomposition as a special case. Now observe 
that 



l<i<i'<k ' ' \<i<i'<k ^ ' 

- ^"Z.Z.^^.i) -'^"Z. -'^fc I n[n-\) ]-"■ 

Therefore, E[i3„] = if and only if a\ =- 0. These results allow us to construct a 
test for (1.2) based on Bn (suitably standardized) as we show in the sequel. 

Firstly using the lexicographically ordered observations, we show that _B„ may 
be reexpressed as 

(2.5) ^""(2) 51 '7«r.^(^r,n), 

^ ' \<r<s<n 

where '\\}{x\^x-2) = [x\ — [i){xi — fi) and 

{n — rii) /{rii — 1), if Yr and Ys arc both observed 

(2.6) rjnrs = \ under the i-th treatment, 
-1, otherwise. 

Thus, it follows that 

(2.7) ^"-^ = 



l<r<s<n 



and 



(2.8) 

1 < r < s < n 



I ■ 



When k 00 and the n^'s are bounded, i.e., the maximum number of observations 
per treatment does not increase with the number of treatments, it follows that 
k = 0{n); hence, by (2.8) we may conclude that M„ = J2i<r<s<n ''Inrs ~ 0{n'^). 

Under Tip : cr^ — 0, the random variables Yi, . . . ,y„ are independent and iden- 
tically distributed and ^j{x,y) is a first-order stationary kernel, centered at 0, con- 
stituting an orthogonal system in the sense adopted in Pinhciro ct al. [19], that 
is, 

(2.9) Vi(n) = m'iYr,Y,)\Yr] = E[iYr - M)(n - fJ-Wr] = a.s. 

(2.10) EmYuY2MY,,Y3)]=0. 

Now, observe that E[i/;2(r^,i;)] ^ YaT[ij(Yr,Y,)] = E[(Yr - fi)'^ {Y, ~ fi)^] = < oo. 
Then, we may show (see Appendix), that under Tig : cr^ = 0, 

(2.11) J„ = -1^^AA(0,1) 

when k oo(=> n ^ oo). The key to the proof is the martingale property exhibited 
by Bn, as pointed in Pinhciro et al. [19] in a different setup. Let lim„^oo Mn/n^ = 
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A, which is finite because M„ = 0{ti?)] assume that the fourth moment of the 
distribution of the random effects is finite and let 5 denote a constant. Then, under 
the sequence of local hypotheses 



Details are presented in the Appendix. Because the mean of this limiting normal 
distribution is positive, we may use J„ as a test statistic for (1.2), rejecting the null 
hypothesis Tio at significance level a when J„ > -Zq, where is the (1 — a) 100% 
percentile of the standard normal distribution. By (2.13), the power of this test 
is directly related to the magnitude of the intraclass correlation coefficient p = 
^b/i'^b + ^e)' more specifically, the power is a monotone increasing function of p, 
as expected. 

3. Simulation results 

We summarize some simulation studies conducted with the objective of evaluating 
the behaviour of the proposed test. First, we examine the efficiency of the J7-test 
in balanced studies under different distributions for the within-treatment errors 
eij and random effects bi. Then, we compare the efficiency of the proposed /7-test 
with that of the usual F-test under various settings. Additional results from the 
simulation studies may be obtained from the authors. 

3.1. Efficiency of the proposed test 

To evaluate the behaviour of the proposed test for small and moderate samples we 
considered 10,000 Monte-Carlo samples obtained under model (1.1) with fj, = 2 
and (Tg = 1, for different numbers of treatments {k — 10,30 and 100) in balanced 
studies. We assumed that ~ A/'(0, 1) and bi ^ x y/afJS, so that Var[6i] = cr^. 
The bctween-treatments variance, a^, was set to (to estimate the size of the test), 
0.2, 0.5 and 1. The empirical power of the test under each setting was evaluated for 
a significance level of a = 0.05. 

The results, displayed in first half of Table 1, suggest that the [/-test tends to be 
liberal when there are few observations for each treatment {rii < 4) . The difference 
between the empirical and nominal size of the test is acceptable when the number 
of treatments is at least 30 and there are 5 or more observations per treatment. 
Furthermore, the power of the test increases with af, rii and k, as expected. 

To illustrate the robustness of the [/-test with respect to heavy-tailed within- 
treatment error distributions, we repeated the simulation study described previ- 
ously assuming that ~ ts x -^3/5 (so that Var[eij] = 1). The results, summa- 
rized in the second half of Table 1, are not considerably different from the previous 
ones; here, however, the empirical size of the test is closer to the nominal size and 
the test also tends to be more powerful than when the within-treatment errors are 
normal. 



(2.12) 




it follows that 



(2.13) 
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Table 1 

Rejection rates (%) for the 5% level U-test in balanced designs with bi ^ tg X 
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3.2. Comparison between the F- and U -tests 

To compare the two tests, we first repeated the simulation study described previ- 
ously, using the same values for /i and in the following situations: 

(i) Balanced study under normality; 

(ii) Moderately unbalanced study under normality; 

(iii) Lightly unbalanced study with heavy-tailed distributions; 

(iv) Balanced study with heavy-tailed asymmetric distributions. 

Firstly, we consider a balanced setting with 5 observations per treatment. The 
random effects and the within-treatment errors were generated according to: bi ^ 
A/'(0, (jf) and Cij ^ Af{Q,l). The results are presented in Table 2 and as in the 
previous studies, they suggest that the [/-test is liberal. However, the difference 
between the empirical and nominal size of the test decreases substantially when 
the number of treatments increases. This occurs mainly when there are few units 
per treatment. On the other hand, the proposed i7-test is more powerful than the 
exact F-test, though we should recognize that this might be influenced by its liberal 
nature. 

We also considered simulation studies to evaluate the effect of imbalance in the 
power of the two tests under investigation. Khuri et al. [IS], Donner and Koval [7] 
and Lee [f ■'3] discuss the performance of F-tests for different degrees of imbalance, 
defined as k = 1/(1-1- -0^) where ip denotes the coefficient of variation associated 
with the sample sizes ni, . . . ,nk- By definition, < k < 1 with k = 1 only for 
balanced studies. Smaller values of k indicate larger degrees of imbalance. 

We selected the number of within-treatment observations from a geometric dis- 
tribution (shifted to the right by 2) with parameter p = 0.15 and support {0, 1, . . .}. 
This corresponds to a moderate degree of imbalance (k = 0.61). The random effects 
bi and the within-treatment errors Cij were respectively generated from A/'(0, cr^) 
and Af{0, 1) distributions. 

In spite of the moderate imbalance, the size of the F-test is closer to the nom- 
inal level than is that of the (more liberal) [/-test, especially when the number of 
treatments is small; however, the latter is more powerful than the F-test in most 
settings. 

We also evaluated the simultaneous effects of a small degree of imbalance and 
heavier-tailed distributions for the random effects and within-treatment errors. In 
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Table 2 

Rejection rates (%) for the 5% level F and U-tests 
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this direction, we selected the number of within-treatment observations with equal 
probability from the set {5, 6, 7, 8, 9, 10}; this corresponds to k = 0.95. In addition, 
we assumed that bi t^.ix ■\/o-^21/41 and eij i4.i x ■\/21/41, so that Var[6i] = cr^ 
and Var[eij] = 1. 

The results, also presented in Table 2, suggest that the size of the [/-test is closer 
to the nominal level than is that of the F-test even in settings with few treatments. 
The bias is larger for the i^-tcst, mainly when there are few treatments. 

Finally, we evaluated the possible effect of asymmetric heavy-tailed distributions 
for the random effects and within-treatment errors on the size and power of the 
two tests. In this direction, wc considered a balanced setting with 5 observations 
per treatment. The random effects and the within-treatment errors were generated 
according to the following distributions: bi ~ {{Yi — E[Yi])/ ■\/Var[Yi]} x erf, and 
Cij ^ {Y2 — E[>2])/-\/Var[y2], where Fi,l2 are independent identically distributed 
random variables with skew t distribution with 4.1 degrees of freedom, location 
parameter 0, dispersion parameter 1 and asymmetry parameter A = 1 with index 
of skcwness equal to 1.77. For details on the skew t distribution, see Azzalini and 
Capitanio [2]. The results are also presented in Table 2 and suggest a conservative 
behaviour for the i^-test. On the other hand, the size of the [/-test is closer to the 
nominal level, even for samples of moderate size. It may also be observed that the 
[/-test is more powerful in all configurations. 

4. Discussion and conclusion 

Although there exists an exact F-test with optimal properties for testing the signif- 
icance of the between-treatments variance component in a one-way random effects 
model with balanced data under normality, we must rely on sub-optimal or approx- 
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imate tests in unbalanced or nonnormal settings. We derived an asymptotic J7-test 
that may be employed with unbalanced data and does not require a specified form 
for the underlying probability distributions. Simulation studies suggest that it has 
reasonable properties even for moderate and small samples. 

In particular, wc may conclude that the F-tcst is more affected by the lack of 
normality of the random effects and within-treatment errors than by imbalance. 
Furthermore, the [/-test seems to be less sensitive to imbalance and to be more 
powerful than the _F-tcst in general. Such conclusions must be viewed with caution, 
given the liberal nature of the [/-test, specially for small sample sizes. To bypass this 
problem, one could consider bootstrap methods to obtain the empirical distribution 
of the statistic J„ under Ho and use the fact that the null hypothesis should be 
rejected for high values of J„. For details on the use of bootstrap to construct 
confidence regions or for hypothesis testing, see Davison and Hinkley [/i] or Lchmann 
and Romano [Ki], Chapter 15, for example. 

In summary, under normality or under moderate degrees of imbalance, the i^-test 
behaves well in terms of power when compared to [/-test; however, in situations 
where the distribution of the random effects and within-treatment errors are non- 
normal, the [/-test is preferable even when the number of treatments is small. In all 
settings the [/-test is more powerful than the _F-test, mainly for small and moderate 
samples. 

The proposed U -test may also be used to test for simultaneous significance of all 
random components in models with more than one random effect, since it behaves as 
(2.11) under Tio. The corresponding power computations, however, require further 
research. 

Brownie and Boos [4] and Boos and Brownie [^] study the behaviour of i^-tests 
and analogous rank tests in one-way random effects models when the number of 
treatments is large. Akritas and Arnold [1] also study the asymptotics of ANOVA 
under similar conditions but their results require more restrictive assumptions than 
those we considered; the proofs, however, are different for balanced and unbalanced 
studies. 

In the context under investigation, the derivation of tests for (1.2) may not fol- 
low the standard procedures since the null hypothesis defines a point (or region) on 
the boundary of the parameter space and this brings in some technical difficulties. 
Asymptotic tests for (1.2) or, more generally, for testing the significance of variance 
components under linear mixed models are available in the literature. Based on 
the ideas of Silvapulle and Silvapulle [35], Verbeke and Molenberghs [38] obtained 
score-type tests under the assumption that the underlying probability distributions 
are normal. Along the same lines, Savalli et al. [20] extended the results to accom- 
modate elliptical underlying distributions. In particular, for the one-way random 
effects models, the corresponding test statistic follows an asymptotic distribution 
given by a 50:50 mixture of Xo and Xi distributions. Tests based on generalized 
likelihood methods (that are asymptotically equivalent to the score-type tests) are 
considered in Self and Liang [22], Stram and Lee [37] and Silvapulle and Sen [36], 
for example. The main disadvantage of such tests is the difficulty in verifying the 
required regularity conditions as shown in Giampaoli and Singer [9] . The derivation 
of the proposed [/-test is not affected by such difficulties and we envisage that it 
may serve as a building block for more general setup as indicated in Nobre [18]. 

Other alternatives have been suggested in the literature. Using a Laplace expan- 
sion of the integrated log-likelihood, Lin [17] obtained a global score test for the 
hypothesis that all variance components are zero in the framework of generalized 
linear models with random effects. Along the same lines, Zhu and Fung [39] ob- 
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tained a global score-type test for the hypothesis that all variance components are 
zero in a semiparametric mixed model, under the assumption of normality only of 
the conditional error vector. It is important to point out, that the tests are ob- 
tained without specifying the distribution of the random effects, though they are 
based on a marginal approach under which the within-trcatmcnt covariances may 
be negative. Hall and Praestgaard [10] consider such a constraint in the framework 
of generalized linear models with random effects, via a projected score test. In prac- 
tice, their results are difficult to apply when the dimension of the vector of random 
effects is large. 



Appendix A: Proof of (2.11) 

Note that under TLq : ct? = 0, 



-1 

n 



l<r<s<7i 



is a sum of (") uncorrelated terms such that E[?/;(yj,, i^s)] = e \a,r[ij;{Yr^Yg 
at < oo. Thus, it follows that 



-2 

n ^ 



(A.l) E[B„] = and Var[B„] = a^Mn 

Then, along the lines adopted in Pinhciro et al. [19], we explore the martingale struc- 
ture of Bn, and apply a martingale Central Limit Theorem as given in Dvorctzky 
[n] or Sen and Singer [33] to obtain the desired asymptotic distribution. Their proof 
requires k fixed and rii oo; we consider rii bounded and let fc ^ oo, although the 
test is also valid when both k oo and — > oo. 
Initially, consider the statistic 

(A.2) r„= Vmj<l>{X^,Xj) 

l<i<j<n 

with (f) satisfying 

(A.3) 0i(Xi) =E[0(Xi,X2)lXi] =0 a.s., 

(A.4) E[0(Xi,X2)0(Xi,X3)] =0, 

(A.S) E[02(Xi,A2)] <oo, 

where Ai , . . . , A„ represents a sequence of independent and identically distributed 
random variables and rjnij denotes weights such that 

(A.6) J2 Vmj^O and ^ '7™, - 

l<i<j<n i<i<j<n 

Defining, Znj = Y^lZi 'nm3 4'iXi, Xj), for j = 2,...n and T„„ = Z)"=2 ^'^^ 2 < 
u < n, it follows that T„„ = T!„. Now, considering the nondecreasing sequence (in 
u) of (T-fields Gnu = cr{Xi, i < u) that corresponds to the cr-fields generated by the 
random vector (Ai, . . . , A„)^, 2 < u < n we may apply the following result due to 
Pinhciro et al. [19]. 
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Theorem A.l. Under conditions (A.2)-(A.5), {r„„,2 < u < n] is a zero mean 
martingale process adapted to the filter {cr„„, 2 < m < n}. 

Letting Dn denote a weakly consistent estimator of = E[(/)^(Xi, Xj)] < oo 
{i < j), it follows from (A.4) and (A.5) that Var[r„] = E[T^] ^ t^M*. Using 
Slutsky's theorem, we obtain 

(A.7) ^^1. 

Now, note that the statistic (2)^11 is of the form (A. 2) with (j){x,y) — (x — fi){y — ^) 
satisfying (A.3)-(A.5). Additionally, observe that = E[4>'^{Xi, Xj)] = < cx). A 
consistent estimator of may be obtained from the following result. 
Lemma A.l. Consider a linear combination of U -statistics of the form Wn = 

T!:=iitu^^T!:=i^^u,. If 

(A.8) max — ^ > 0, 

it follows that W^-^a^. 

Proof. Condition (A.8) is valid when the n.^'s are fixed and /c ^ 00. It also holds 
when n,; — > 00, as for example, in a balanced study. 

Noting that E[W,i] = ct^, and using properties of ?7-statistics, it is possible to 
show that Var[[/i] = E[e|j]/ni — (n^ — 3)CTg/{(ni — l)ni}. Then, a direct application 

of Chebyshev's inequality leads to Wn—^cr'e- Since is a continuous function in 

the support of Wn, we conclude that Dn = W^j-^crf. 

In the proof, E[efj] < 00 is used to guarantee the existence of Va.r[Ui\. In some 
situations, as in the balanced case, this requirement may be relaxed and the result 
follows from Khintchine's weak law of large numbers. □ 

We now focus on the following lemma. 
Lemma A. 2. Consider the weights defined in (2.6). When k ^ 00 it follows that 

(A.9) vlruVlju/M^^^O, and ^ ?4/A/,^ ^ 0. 



Proof. When k ^ 00 and the n^'s are bounded, it follows that k ~ 0{n). Then by 
(2.8), we obtain A/„ = 'l2i<r<s<n''lnrs = 0{n^). Now, the result follows from 

^ntj = Oin^) and ^ vltuVlju ^ 0{n^). 

^^i^j^n l<.i^j <u<.n 

□ 

The proof of the main result consists in verifying the two regularity conditions 
of Theorem 3.3.7 in (Sen and Singer [:>.!], page 120) and it is similar to the proof 
presented in Pinhciro ct al. [!!)]. 

Theorem A.2. Suppose that E[0'*(Ai, A2)] < 00, (A.3)-(A.5) hold and, for n 
00, the weights rjnij satisfy 

(A.IO) Y vLkVl,k/{M*)'^0 and ^ vLj/{M*)'^0. 

l<.i^j <k<.n l^^T^J^""- 
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Then, for r„ in (A. 2), we have that 

(A.ii) (A/,:A^)-'/'rn^AA(o,i). 



Proof. For the weights defined in (2.6), we have M* = Then (2.11) follows 
from Lemmas A.l and A. 2. □ 

Appendix B: Proof of (2.13) 

First observe that 

i<j<j'<«. ^ i<j<i'<n. ^ 

fc 

and W,i = 'y^^{ni/n)Ui, are functions only of the vector of within-treatment er- 

i=l 

rors e. Consequently, their distributions under Tip and 7ii„ are the same. On the 
other hand, we have 

i?„ = V -^^^A2Uu'-U..-Ue} 
^-^ n(n — 1) 

l<i<i'<k ^ ' 

n^n^' I 2 y> {Y,j -Yyyf _u 

<i<i'<k ^ ' I j = lj' = l 



1< 



where 



^^4^^ n(n-l) n,n,,^^ 2 f ' 



a. = E (^.-^•^o^ + ^^^^EE(e.-e.,o 

71(71— 1) 71i7li' ^-^ ^-^ 

l<i<i'<k I, j=ij' = i 

Also, under Hq, -B„ = B^. Therefore, 
m n r - (2)^" _ (2)^° , (2)'^" _ . , ^ 

where J„o = ('20i?°/('^e ^^r^ )^AA(0, 1) (by 2.11) and Q„ = (;')C„/(a2A/y2). 
The term C„ can be decomposed as C„ = C„i + C„2 where 

C„i = , ^ niU^' {bi - h')'^ 

n(n — 1) 

^ ' l<i<i'<k 

and 

" 7^(71 - 1\ 5Z (^^-^'OEE^^'^^^'j')- 

^ ' l<i<i'<fe 3 = 1 j' = l 
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Since E[C„2] = 0, we have 



nc'ni] 



n(n — 1) 

^ ' l<i<i'<k 



n{n- 1) ^ 



<i<i'<k 



71(71 — 1) 



The following result provides the asymptotic behaviour of Q„ under the assumption 
of the existence of the fourth moment of the random effects. 

Lemma B.l. Consider model (1.1) under Hin = 0-5= S'^/y/n. Assuming that 
E[bf] < 00 and letting lim„^oo Mn/n'^ = X, we have 



2^/A' 



Proof. Given that E,[bf] < 00, then by contiguity of the sequence of local alternative 
hypotheses Hin ■ crl = (^^/v^i '^^ have that E[6f] = o(l) under 7ii„. Defining 
fo*-' = {bi — bjY, recalling that bi, . . . ,bk are independent and identically distributed 
and that the 71,; 's are bounded, we have 



Var[C„ 



1 



71^(71 — 1)^ 



E 

l<i<i'<k 



n^nlYarlb'' 



l<i<t<k 
l<i<w<k 



7i^7ijn^Cov[6'*, 6'' 



max{Var[6i2l,Cov[6i2,6i3i} 2 



n'^{n - 1)^ 



0(1)0(71-1) =o(7i-i), 



l<i<i'<k 



l<i<t<k 
l<i<w<k 
Mm 



since Var[&;^^] and Cov[6;!,^, 6;^^] arc of the same order as E[fe^], namely, o(l). Sim- 
ilarly, it is possible to show that Var[C„2] = 0(71-^). Here, observing that (2)/ 
y/ Mn = 0{y/n), we have 



lim 



= lim 



and since 0(7i)Var[C„i] = o(l) 0, when 71 00, it follows that [^)Cni/ 
VWi^ Since E[C„2] = and Var[C„2] = 0(77-1), j^^^^g ^-^^^ 

{t)Cn2/VM;:^^- Therefore, 



C„/v^ = Q (^"1 + C„2)/v^^^V{2^/A}. 



□ 



Now, recalling (B.l), Lemma B.l observing that Wn — ^ol, and using Slutsky's 
theorem, we obtain 
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implying that, under 7ii„ : cr^ = j ^fn the center of the normal is shifted to the 
right by 5V(2cr2VA). 
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